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SUMMARY 


A wind-tunnel investigation was made to determine beat-transfer 
distributions on three steel sphere -ellipsoid bodies with surface rough- 
nesses of 5, 100, and 200 microinches. Tests were conducted in the Langley 
9- by 6 -foot thermal structures tunnel at a Mach number of 3 . 0 , free -stream 
Reynolds numbers (based on model spherical diameter) of 4.25 X 10^ and 
2.76 X 10°, and at a stagnation temperature of 650 ° F. Pressure distribu- 
tions were obtained also on a fourth model. The results indicated that 
the combination of surface roughness and boundary-layer cooling tended to 
promote early transition and nullify the advantages attributable to the 
blunt shape of the model for reducing local temperatures. Good correla- 
tion between experimental heating rates and those calculated from laminar 
theory was achieved up to the start of boundary-layer transition. The 
correlation also was good with the values predicted by turbulent theory 
for surface stations downstream from the 45° station. 


INTRODUCTION 


The present investigation was conducted to determine the effects 
of surface roughness on the heat -transfer distribution over the plastic 
sphere -ellipsoid models of reference 1. In that investigation, wind- 
tunnel air stream contamination roughened the model surfaces to a degree 
which was believed to have affected the location of boundary-layer tran- 
sition and, hence, the region of maximum aerodynamic heating input. 
Inasmuch as instrumentation limitations imposed on the plastic models 
by structural considerations precluded direct measurement of the heat- 
transfer distributions, steel models having different surface roughnesses 
were constructed and tested. 
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Four models were tested simultaneously in the Langley 9- hy 6 -foot 
thermal structures tunnel at a Mach number of 3-0, at free-stream Reynolds 
numbers (based on model spherical diameter) of 4.25 X 10^ and 2.76 X 10 , 
and at a stagnation temperature of 65 O 0 F. Three of the models were 
instrumented with thermocouples in order to obtain heat-transfer data, and 
the fourth model was used to obtain pressure distributions. The pressure- 
distribution model and one heat-transfer model had a surface finish of 
5 microinches. The other heat-transfer models were finished to average 
-surface roughness values of 100 and 200 microinches, respectively, in 
order to cover the range of surface roughness conditions estimated for 
the plastic models of reference 1. A discussion and analysis of the 
results are presented herein. 
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pressure coefficient 
specific heat of model wall material 
specific heat of air at constant pressure 
model spherical diameter 
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thermal conductivity for air 
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Reynolds number whose characteristic length is D, 
local Reynolds number, 


hD 

k 

c pP 

IT 

_h 

PV C P 



(4 H N^O\IA 


3 


r local radius perpendicular to model longitudinal axis 

s surface distance from stagnation point 

T temperature 

t time 

' v velocity 

w specific weight of model wall material 

x axial distance from stagnation point 

0 velocity gradient along the surface, 

c)s 

9 angular distance around models measured from the stagnation 

point in spherical polar coordinate system 

H absolute viscosity of air 

p mass density of air 

t model wall thickness 

Subscripts: 

D parameter based on model spherical diameter as characteristic 

length 

1 local free-stream conditions evaluated at outer edge of 

boundary layer 

o initial conditions 

r recovery, surface conditions 

s parameter based on surface distance as characteristic length 

t stagnation conditions 

w outside wall of model 

a free-stream conditions evaluated immediately behind bow 

shock wave 

oo free-stream conditions ahead of shock wave 
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SPECIMENS, APPARATUS, AND TESTS 


Models 

Principal model dimensions are shown in figure 1. The models were 
machined to the same size and profile as the plastic sphere-ellipsoid 
models of reference 1. Integral with the models was a conical skirt 
with a flare angle of approximately 5° which extended to a base diameter 
of 5-75 inches. Model wall thickness was 0.030 inch for the heat-transfer 
models and approximately 0.12 inch for the pressure model. Heat-transfer 
model wall -thickness variations were held within ± 0.001 inch. The outer 
surfaces of two of the heat -transfer models were sandblasted to average 
roughnesses of approximately 100 and 200 microinches, respectively. The 
outer surfaces of the third heat-transfer model and the pressure model 
were polished to a finish of approximately 5 microinches. Average rough- 
ness heights of the model surfaces were measured with a Physicists Research 
Company Profilometer. In order to prevent large changes in model surface 
roughness during the tests, the models were machined from a hard material 
(SAE 1020 carbon steel) . Precautions were also taken to minimize the 
amount of foreign particles in the airstream by the installation of a fine 
mesh screen at the upstream end of the tunnel nozzle entrance cone. The 
nozzle and test-section areas were also vacuum cleaned prior to each test 
to remove dust and sand that might have been sucked into the tunnel during 
the shutdown period of the previous test. 

Instrumentation for the heat-transfer models consisted of l4 
No. 30 gage iron constantan thermocouples spot-welded to the inner sur- 
face of each model. The thermocouples were staggered along two meridians 
l80° apart emanating from the stagnation point. Thermocouple surface 
station locations are tabulated in figure 1. The pres sure -distribution 
model contained 18 pressure orifices 0.070 inch in diameter. Eight 
orifices were arranged in a spiral of one revolution around the model 
and were located at the surface stations tabulated in figure 1. At 
stations 1.00 and 2.27 and on the conical skirt, there were four orifices 
equally spaced around the circumference. Model pressures and temperatures 
were recorded by means of a high-speed digital magnetic tape recording 
system. 


Test Facility 

The tests were conducted in the Langley 9- by 6-foot thermal 
structures tunnel, an intermittent supersonic blowdown facility exhausting 
to the atmosphere. Air for the facility is dried to a dewpoint of -40° F 
and stored in a large tank field at a pressure of 450 pounds per square 
inch absolute. Four quick-acting rotary valves automatically control the 
air from the storage field to provide preset stagnation pressures from a 
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minimum starting pressure of 120 pounds per square inch absolute to a 
maximum of 200 pounds per square inch absolute. Stagnation temperatures 
up to approximately 660° F are obtained by passing the air through a 
stainless- steel heat exchanger which is preheated by an oil-fired indirect 
air heater. The air then passes through a two-dimensional supersonic 
nozzle designed for a test section Mach number of 3*0. 

The test section is 6 feet high, 8 feet 9 inches wide, and 10 feet 
long. Calibration tests have shown that the maximum deviation in test 
section Mach number is less than 1.0 percent. Flow conditions and running 
times available in the test section are given in the following table: 


Mach number 

6.6 > 

3.0 

Reynolds number per foot 

< 10 6 to 19 x 10 6 

Dynamic pressure, lb/sq ft 

3,000 to 5,000 i 

Stagnation temperature, °F 

300 to 660 

300 to 66o 

Maximum running times for constant Tt, sec . . 

23 

28 

15 

18 

Maximum running times for constant p t , sec . . 

4o 

1+8 

18 

22 


The lower stagnation temperature limit of 300° F is imposed to 
avoid the possibility of condensation in the test section. The running 
times listed are the maximum times for which either the stagnation pres- 
sure or stagnation temperature can be held constant. Running times for 
which stagnation temperature is held constant are based on a temperature 
drop of not more than 10° F. All times shown in the table are exclusive 
of the starting and shutdown periods of 2 seconds and 7 seconds, respec- 
tively. During these periods, the flow separates from the nozzle walls 
and subjects the test specimens to unsteady loads considerably in excess 
of the aerodynamic loads applied under test conditions. 


Model Installation 

All four models were tested simultaneously. Each model was attached 
to a cylindrical mounting bucket at the end of a short strut as shown in 
figure 2. The heat -transfer models were mounted from the upper test- 
section wall with the thermocouples oriented in a vertical plane. The 
pressure model was mounted off a side wall. All models were located at 
the same tunnel longitudinal station at zero angle of attack and yaw 
with respect to the tunnel longitudinal center line. 
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Test Conditions 

Two tests were performed. Test conditions were selected to cover 
the limits of the range of conditions imposed during the tests of ref- 
erence 1. Free-stream Reynolds numbers (based on model spherical diam- 
eter) were 4.25 X 1C>6 and 2.7 6 X 10^. Tunnel stagnation pressures were 
202 and 152 lb/sq in. abs, and the stagnation temperature was 650° F. 

The duration of test conditions for the high Reynolds number test was 
18.5 seconds, whereas the lower Reynolds number permitted extending 
this period to 26.2 seconds. These values are exclusive of the starting 
and shutdown periods of 2 seconds and 7 seconds, respectively. 


RESULTS AMD DISCUSSION 


The condition of the model surfaces after the tests showed that the 
airstream contamination had largely been eliminated. Profilometer readings 
taken on the heat-transfer and pressure model with an initial surface 
roughness of 5 microinches showed that the surface finish changed by only 
20 microinches. 


Model Pressures 

Pressures measured along the sphere-ellipsoid surface and conical 
skirt at the two free-stream Reynolds numbers of this investigation are 
shown in figure 5 as the ratio of local to stagnation -point pressure. 
The measured pressures are compared with the pressures calculated from 
modified Newtonian theory given by the relation 


Ct, = C. 


p,a 


cos^9 


As shown in figure 5 > the measured pressures are closely approximated 
by the theoretical values over most of the surface. Consequently, 
isentropic flow relations and the modified Newtonian pressure distribu- 
tion were used in calculating the local flow conditions required in the 
analysis of the heat -transfer data discussed in the appendix. 


Temperature Distributions 

Temperature data were obtained from all three heat-transfer models 
during the high Reynolds number test, but extensive instrumentation 
system failures occurred in both the smooth model and the model with a 
surface roughness of 100 microinches during the low Reynolds number 
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test. Consequently, only the data obtained from the model with a sur- 
face roughness of 200 microinches are reported for the low Reynolds num- 
ber test. 

Typical surface -temperature distributions obtained from each model 
at various times during the high Reynolds number test are shown in fig- 
ure 4. The sharp rise in temperature noted downstream from the stagna- 
tion point in figure 4 is characteristic of boundary-layer transition 
to turbulent flow. As the test progressed, the peak in the curves 
gradually diminished to the relatively flat distribution shown for 
t = 18 seconds. In general, equilibrium temperatures were achieved just 
before the end of test conditions. 


The Effect of Surface Roughness on 
Heat-Transfer Distribution 

Experimental heat -transfer results obtained early in the tests from 
each model are plotted in figure 5 as the distribution along the surface 

of the laminar heat -transfer parameter, ^Rp^. The airflow prop- 

erties used were the undisturbed free-stream conditions, and model 
spherical diameter was taken as the characteristic length. Also shown 
in figure 5 are theoretical curves for both laminar and turbulent boundary 
layers. The curves were calculated from the modifications of Sibulkin's 
and Falkner's theories given In references 2 and 3- The procedures used 
in determining the experimental heat-transfer coefficients and theoretical 
values of the parameter are discussed in the appendix. Errors in the 
experimental heat-transfer coefficients introduced by conduction effects 
were evaluated at surface stations where these effects were considered 
to be greatest by using the method given in reference 3» Since the maxi- 
mum change in heat-transfer coefficient at these stations varied from 
approximately 6 to 8 percent and the average change at other stations was 
less than 4 percent, the effects of conduction were neglected. Heat 
losses due to radiation were also evaluated and considered small enough 
to neglect. 

The experimental heat-transfer distributions shown in figure 5 
indicate that all three models experienced transition to turbulent flow 
as evidenced by the abrupt departure from the trend of the laminar theory 
curve. The most evident effect of surface roughness is the familiar 
upstream shift in the location of boundary-layer transition. For these 
models, the start of transition moved from the vicinity of the sphere- 
ellipsoid junction (a station 45° from the stagnation point) on the 
smooth model to a position just downstream from the stagnation point 
on the roughest model. The indicated shift in the location of the start 
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of transition between the models with surface roughnesses of 5 and 
100 micro inches was greater than the shift shown between the models with 
surface roughnesses of 100 and 200 microinches. Accompanying this upstream 
shift in the location of transition due to roughness is an increase in peak 
heating rates. For the model with a surface roughness of 200 microinches, 
these heating rates attained values approximately twice the stagnation- 
point value. 

The good agreement between the experimental values and the values 
obtained from laminar theory up to the start of transition, as seen in 
figure 5, indicates the adequacy of the theory for predicting laminar 
heat-transfer distributions for the bodies tested in this investigation. 
Also, the agreement between experimental values and the values given by 
turbulent theory is good at surface stations downstream from the sphere - 
ellipsoid junction. 

In the present investigation, stagnation -point values of the 
parameter, Ng^ ^R-^ ~ , were not obtained for the two roughest models 

(see fig. 5) , and a comparison cannot be made to ascertain differences 
in stagnation-point heating due to roughness. However, a direct compari- 
son of the values of the heat -transfer parameter on the three models 
can be made in figure 5(a) at surface station, ^ = 0.0 56, just down- 
stream from the stagnation point. Although the value of the param- 
eter, oo, is highest for the roughest model, this may indi- 

cate only that the flow is not laminar at that station. Reference 
to figures 5(a) and 5(b) shows an increase in the heat-transfer parameter 
with increase in free -stream Reynolds number for the roughest model 

(200 microinches surface roughness) at surface station, ^ = O.O56. 

In reference h, an increase in the laminar heating rate with increase 
in free-stream Reynolds number was observed on a roughened hemisphere- 
cone -cylinder. For the present investigation, however, there is no 
assurance that an increase in the laminar heating rate is represented 
for the test made at the higher Reynolds number. 


Boundary-Layer Transition Reversal 

In determining experimental local heat-transfer coefficients, two 
values were obtained for each model at surface stations in the vicinity 
of initial peak heating rates shown in figure 5- For these stations, a 
decrease in heat -transfer coefficient occurred when the model wall tem- 
perature attained a certain value. The wall temperature at which the 
reduction in heat -transfer coefficient took place increased as the sur- 
face distance measured from the stagnation point increased. This effect 
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is illustrated in figure 6 , "where the wall-to-free-stream temperature 
ratio for reduction in heating rate is plotted against free-stream 
Reynolds number with surface distance, s, as characteristic length. 

As shown in figure 6, the wall-to-free-stream temperature ratio varied 
from approximately 2.33 to 2.67, depending on the surface roughness and 
surface distance. For the model with a surface roughness of 200 micro- 
inches, the wall temperature ratio for reduction in heating rate was 
found independent of the free-stream Reynolds number based on model 
spherical diameter and also appeared to approach a constant value. 

The change in heat-transfer distributions due to the lower heating 
rates is shown in figure 7? where the values of the param- 
eter, Ng-k j00 \jRp j0 o, obtained later in the tests are shown as the closed 
symbols. For the smooth model (5 micro inches surface roughness) at 
Rp,oo = 4.25 X 10^ (fig. 7(c) ) and the roughest model (200 microinches 
surface roughness) at Rp „ = 2.76 X 10^ (fig. 7(d)), the reduction in 
heat-transfer coefficient was sufficient to conclude that the flow had 
changed from turbulent to laminar. Thus, the transition point moved 
downstream after an increase in model wall temperature. A comparison 
of the conditions at which the downstream movement of transition first 
occurred on the smooth model tested at the high Reynolds number (see 

fig. 7(c), surface station ^ = 0.333) with those of the roughest model 
tested at the lower Reynolds number (see fig, 7(d), surface sta- 
tion H = O.167) showed that the local pressures and model wall tem- 
peratures at these two surface stations were nearly the same. 

A downstream shift in the location of transition on smooth blunt 
bodies with increase in surface temperature was also observed during 
tests reported in references 4 to 6. This behavior, known as transition 
reversal, has been the subject of considerable interest (see, for example, 
refs. 4 to 9 ) y and its cause is not yet thoroughly understood. The chief 
point of controversy is centered on whether the apparent boundary-layer 
destabilizing effect of extreme surface cooling is of a nonroughness 
origin. (See, for example, ref. 8.) In reference 6, it was concluded, 
on the basis of calculations of the minimum roughness height required to 
excite transition, that roughness was not the dominant factor in causing 
the initially turbulent flow during the hemisphere tests of that investi- 
gation. On the other hand, in reference 7 it was shown experimentally 
that transition can result directly from the effects of roughness in a 
cooled boundary layer. For the present investigation, it is clear that 
the combination of surface roughness and boundary-layer cooling tended 
to promote early transition to turbulent flow and nullify the advantages 
attributed to the bluntness of the model shape in reducing local 
temperatures. 
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CONCLUSIONS 


An investigation to determine the heat-transfer distribution on a 
sphere -ell ip so id body for three surface roughness conditions was con- 
ducted in the Langley 9- by 6-foot thermal structures tunnel. Three 
thermocouple instrumented, steel models were constructed with surface 
finishes of 5, 100, and 200 microinches, respectively. The models were 
tested simultaneously at Mach number 3«0, free -stream Reynolds numbers 
(based on model spherical diameter) of 4.25 X 10^ and 2.76 x 10^, and 
a stagnation temperature of approximately 65 O 0 F. Local pressures were 
measured on a fourth model. The test results indicated the following: 

1. As surface roughness increased, the location of the start of 
boundary-layer transition shifted from a station 45° from the stagnation 
point on the smooth model (5 micro inches surface roughness) to a position 
just downstream from the stagnation point on the roughest model 

(200 microinches surface roughness). A greater difference in the location 
of transition was observed between the smooth model (5 micro inches surface 
roughness) and the model with a surface roughness of 100 microinches than 
between the two roughest models (100 and 200 micro inches surface roughness) 
The differences in location of the start of transition due to roughness 
were accompanied by an increase in peak heating rates which attained values 
approximately twice the stagnation-point value for the roughest model. 

2. As the model surface temperature reached certain values, a 
decrease in heat-transfer coefficient occurred on all models at surface 
stations in the vicinity of initial peak heating rates. The wall-to- 
free-stream temperature ratio at which the reduction in heating rate 
occurred increased as the surface distance from the stagnation point 
increased and varied from approximately 2.33 to 2 . 67 . The reduction in 
heat-transfer coefficient was sufficient to conclude that the location 
of boundary-layer transition had moved downstream on the smooth model 
for the high Reynolds number test and on the roughest model for the low 
Reynolds number test. 

3 . In general, the measured pressure distributions over the model 
were in good agreement with those predicted by modified Newtonian theory. 

4. Theoretical heat-transfer distributions for laminar flow using 
modifications of Sibulkin's theory were in good agreement with experi- 
mental distributions up to the start of transition. Also, theoretical 
heat -transfer values for turbulent flow using a modification of Falkner's 
theory were in good agreement with experimental values at surface stations 
downstream from the sphere -ellipsoid junction. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Field, Va. , April 13, 1961 . 
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APPENDIX 

DATA REDUCTION AND THEORIES 


Experimental Heat-Transfer Coefficients 

Local heat-transfer coefficients were determined from the slope 
of the straight-line variation of the measured temperature plotted 
against its time rate of change in the heat -balance equation 


c m VT 


dT w 

5T “ 


(i) 


where h and T r are assumed constant for a given set of test condi- 
tions. In equation (l) T r is the intercept on the T w axis when T w 


is plotted against 


dt w 
dt ' 


In order to ensure against using data influ- 


enced by irregularities in test conditions, no data were considered in 
the analysis at times below 2.5 seconds. At some surface stations, two 


straight-line variations of Tw with — — were obtained as shown in 

0.0 

figure 8. Using the slopes of each straight line, two values of h 
were computed. This indicated the existence of the transition reversal 
phenomenon in these tests as shown in figure 7- The value of h obtained 
in this manner was checked by computing temperature variations with time 
at each station from the solution to the differential equation (l) given 
by the relation 


Tw ~ T 0 
T r - T 0 


= 1 


ht 

” c m WT 
e 


( 2 ) 


The difference between computed and measured temperatures was always 
within 2 percent of (T r - T 0 ) . Physical constants used in the data 
reduction were c m = 0.129 Btu/(lb)(°F) and w = ^90 Ib/cu ft. 


Heat-Transfer Parameter 

Stagnation point . - The curves labelled "laminar theory" in figures 5 
and 7 were computed from modifications of Sibulkin's theory given in 
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reference 2 for a hemisphere -cylinder. A stagnation-point solution was 
presented in reference 2 in terms of the Nusselt number and Reynolds 
number based on spherical diameter and physical properties of the air 
evaluated immediately behind the normal shock. The theoretical heat- 
transfer parameter obtained is a function of Prandtl number and velocity 
gradient along the surface as follows: 


%u,D,cr 


= 0 . 7 63 ( 


N- 


Pr,a 


2/5^3Dj!/2 


(3) 


By expressing this solution as the product of the Stanton number and 
the square root of the Reynolds number, equation (3) becomes 


N St,a\j R D,<X 




■ N Pr,tf>PD,ff 


(4) 


Transformation to undisturbed free-stream conditions required multiplying 
both sides of equation (4) by the product. 



p a Y a^a 

Pco^co^oo 


1/2 


Equation (4) can then be written as 

1 / 2 / c P ,a \ f p o Y a^c 

V C P^A P “ Voo4 °o 

( 5 ) 

Using isentropic-flow relations and a value of ^ = 4.97 determined 

from modified Newtonian theory, a stagnation -point value of 
%t,oo-j%i,oo =3*44 is obtained. This value is in agreement with that 

computed at the stagnation point of a hemisphere-cone-cylinder in 
reference 4. 


N St,°o\J R D,c 


PooVoob 


J 
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Laminar flow . - In reference 2, a relation is given for the local 
Stanton number which should he valid for stations up to from the 

stagnation point or £ = 0.393* The relation, also based on a. modifica- 
tion of Sibulkin's theory, uses physical properties of the air evaluated 
at the outer edge of the boundary layer and s as characteristic length 
as follows: 


- °-T63( H P r ,j )' 3/5 {&) 

Transformation to undisturbed free -stream conditions and spherical diam- 
eter as characteristic length yields 


N. 


St 


,°°\ %),00 


P V C 

^00 00 PjOO 


p v D 
K 00 00 


0.763(N Pr , z ) 



(7) 


The laminar theory curves shown in figures 5 and 6 were calculated using 
equations (5) and (7). Local flow conditions were evaluated from 
isentropic-flow relations and the modified Newtonian pressure 
distributions. 


Turbulent flow . - The theoretical heat-transfer distributions for 
turbulent flow shown in figures 5 and 6 were obtained from a modifica- 
tion of Falkner's theory given in reference 3 for a sphere. The rela- 
tion is given as the product of the local Stanton number and \[r s as 
follows: 


%t,Z\P*s = 


h 


P v c, 
w l p ,1 




p w v Z £ 


p. 


w 


0.0157 



5/i4 


(8) 


In this expression, the physical properties of the air are evaluated at 
the wall with s as characteristic length. Transformation to D as 
characteristic length and to undisturbed free-stream conditions yields 
the equation 


%t ,oo^%), 00 


0.0157 





-1/7 


(9) 



For the turbulent theory curves shown In figures 5 and 6, the local 
flow conditions were determined from isentropic-flow relations and the 
modified Newtonian pressure distribution. 
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Figure - Surface temperature distributions obtained at various times 
from three sphere -ellipsoid models finished to different surface 
roughness heightsj times shown are based on time zero equal to 
0.5-second test time; Mach number 5>0; T + = 650° F. 














Figure 6.- Variation of model wall temperature for reduction in heating 
rate with free-stream Reynolds number based on surface distance from 
the stagnation point on sphere -ellipsoid bodies for different surface 
roughness conditions; spherical radius, 2.25 inches. 
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Figure 7- - The effect of increasing wall temperature on the local heat 
transfer to sphere-ellipsoid bodies with different surface roughness 
conditions. 











Figure 8.- Method for determining experimental heat-transfer coefficients 
from thermocouple dataj illustration of transition from turbulent to 
laminar flow of the smooth model at s = 2.273 inches. 
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